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F, Erdogan and M. Bakioglu' 


ABSTRACT: The tracture problem in panels which consist of periodically arranged 

load-carrying and buffer strips of different materials is considered. The main 
nmphasis in the study is placed on the problem of a crack terminating at and 
rirossing the interfaces and on the stress-free end problem. The problem is for- 
mulated in terms of a system of singular integral equations and numerical solu- 
tions are obtained for certain material combinations. With the study of possible 
crack propagation and delamination in mind, certain stress intensity factors are 
defined and calculated, A main result of the study is that when the crack touches 
or intersects a bimaterial interface, the stress state has no longer the standard 
square root singularity and, to study further propagation of the crack, the con- 
ventional fracture models need to be modified or new models need to be developed. 


KEY WORDS: Con^posite materials, fracture, crack propagation, delamination, stress 

singularity, bonded layers, stress-free end. 


In structural design with high strength composite sheet materials the use of 
relatively low stiffness and high toughness buffer strips oriented parallel to 
the main load-carrying laminates has been under investigation for some time [1,2]. 
The practical objective of this design procedure is to improve the fatigue crack 
propagation and arrest characteristics of the structure. The fracture process may 
start as a fatigue crack initiation at a local imperfection in the load-carrying 
laminate. Under repeated loading it may be possible for the propagating crack to 
reach the interface and to enter the adjacent buffer strips. Assuming that the 
fatigue and fracture characteristics of both materials are known, for studies re- 
lating to fatigue life and structural integrity, it would be necessary to have a 
reliable analysis of the problem. Particularly in fatigue crack growth studies, 
it would be very useful to have a technique for calculating the corresponding 
stress intensity factors. 

After the crack propagating in the first strip touches the interface, further 
fracture propagation would take place along the plane for which the ratio of the 
local load factor (or crack driving force) representing the intensity of the 
applied loads and the geometry to the corresponding strength parameter of the two 
materials or the interface bond is maximum. In this rase, sines the stress state 
around the crack tip does not have the standet..d square root singularity, the con- 
ventional fatigue and fractxure models cannot be used to study the problem. To 


'Professor 

^Lecturer, 

Istanbul, 


of Mechanics, Lehigh University, Bethlehem, PA 18015. 
Department of Civil Engineering, Istanbul Technical University, 
Turkey. 


- 1 - 


dc» this, in addition to having the correct analysis of the problem it may be 
necussary to make certain modifications in the existing fatigue and fracture 
models. 

In the actual structures, particularly in those which consist of fiber re- 
inforced laminates, generally the strips are anisotropic [l]. However, in this 
study, largely for reasons of analytical expediency, both materials are assumed 
to be Isotropic and linearly elastic. The basic formulation of the problem of 
periodically arranged bonded dissimilar strips and its solution for cracks grow- 
ing in one set of strips and terminating at the interfaces were given in [2] 
(Figure la) . In this paper the anphasis is on the propagation of the crack in 
the second medium after crossing the interface and on the "stress-free end prob- 
lem" of bonded layered materials (Figures la and lb) . 

Fo^aCcuUon 0^ th& PA.oblem 

The general problem under consideration is described in Figure 1. it is a 
two-dimensional elastostatics problem for a composite medium which consists of 
periodically arranged two sets of bonded dissimilaj* strips having collinear 
cracks as shown in Figures la and lb. The medium juay be under symmetric uniform 
normal tractions in x, y, and (ir- plane strain case) z directions, or under in- 
ternal stresses due to the difference between curing and operating temperatures. 
Under any combination of these loads the problem giving the homogeneous stress 
state in the composite medium in the absence of cracks may be solved and the 
stress components on y=0 plane may be obtained without any difficulty. Let 
these stresses be 


= Pi '-I 


00 

ou 




( 1 ) 


The solution of the original crack problem is obtained by adding to the homo- 
geneous solution the solution found from a perturbation problem in which the 
self-equilibrating crack surface tractions are the only external loads. From 
the view point of fracture the important problem is the perturbation problem 
which may be solved by considering the standard field equations for the strips 
1 and 2 under the following boundary and continuity conditions: 


Uj^(0,y) = u^(0,y) , Vj^(0,y) v^O.y) , (0<y<«>) , 
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Ui(-hj^.y) - 0 , ^ ' (0<y««) , 

U2(h2,y) » 0 , °2xy^^2'^^ “ ' 

0lxy(x,0) = 0 , (-h^<x<0) , «^2xy^’‘'°^ “ ° ' ( 0 <x<h 2 ) , 

Oj^yy(x,0) *= “P]^<x) , (-h^<x<a-hj^) , Vj^(x,0) » 0 , (a-hj^<x<0) , 

O2yy(’*'0) “ -Pjfx) t (h2-b<x<h2~c) , V2(x,0) « 0 , (0<x<h2”bf h2”C<x<h2) » 
°kyy^’^'“^ “ ^ ' ^kxy^*'“^ “ ° ' (2) 


where, because of symmetry, only one quarter of a strip from each set is con- 
sidered, Uj^, Vj^, the dibplacement and stress components in the usual 
notation (k==l,2; i,j®x,y), and the dimensions h^^, h2, a, b, and c are shown 
in Figure 1. Defining 


xi = X + hi , X2 = X - h2 , yi = y - y. 

GilXiJ = Ji|“ , U = i,2) 


(3) 

(4) 


rtud uain>j the t«.c;hnique described tn [2], the perturbation problem may be re- 
duced tr, the following system of singular integral equations: 
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^Gi<t)dt4 
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kii (Xi,t)Gi(t)dt 


iTa+<i) 

k, „(xwt)G-(t)dt T~ p, (x -h.) , ~a<x<a , 
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■ e -2 <=« 2-“2 


•)G.(t)dt + k^, (x.,t)G, (t)dt 


21 ' 2 ' 


+ I k22(X2,t)G2(t)dt 
c 


it(1+<2) 

“Hu 


P2(X2+h2> , b<X2<b , 


(5) 


where ^i, iCi, (i=l,2) are the elastic constant.s with <i= (3-Vi)/ (1+Vi) for 
generalized plane stress, < 1 « 3-4Vi for plane strain, Vi being one Poisson's 
ratio. The kernels kij(Xi,t), (i,j = l,2) are given by 


k.i(x.,t) = 


Kiitx.,t,s)e“®^^l‘^^ds , 


ki2(Xi,t) = I [Ki2(Xi,t,s)e - Ki^fXij-t.sje ®'‘^2^^^]ds (i = l,2) 


( 6 ) 

where the functions are defined in the Appendix B of [ 2 ] and will not be 
repeated in this paper. The index of the integral equations (5) is +1 [ 3 ]. 
Therefore the general solution of (5) will contain two arbitrary constants [ 3 ] 
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which may be determined by using the following single-valuedness conditions: 


-a 


Gj^(t)dt = 0 


G^(t)dt = 0 


( 7 ) 


For cracks imbedded in the strips, i.e., for a<hj^ and 0<c<b<h2, the solu- 
tion of (5) is of the form 

= gj^(t) (a^-t^) , -a<t<a , 

G 2 (t) = g^ (t) [ (b-t) (t-c) ] ** , c<t<b , (8) 


where g^^ and g^ are bounded in their respective closed domains. The unknown 
functions g^^ and g^ may be determined from (5) and (7) in a straightforward 
way by using the technique described in [ 4 ]. For this case the stress intensity 
factors at the crack tips are defined and are expressed in terms of G^^ and G 2 
as follows: 




"a • “5 [2(x^-a)]'’ojyy(Xj,0l . [2(a-Xj)]\(x^> , 


k = lim [2(x -b) 
t> X 2 -^b 2 2yy 2 


1+Kj^ Xj^.^ 

4P 


lim [2(b-x )]*’g (x ) , 


4P. 


k_ = lim [2(c-x.,)] a.,._ (x.,,0) = -z-— lim [ 2 (x_-c) ]^G_ (x_) . 


X2->c 


2' ■' ^2yy' 2 


l+<2 >'2^c 2 ''2'"2' 


(9) 


The case of a crack touching the interface, i.e., a = hj^, was treated in [ 2 ] 
where it was shown that 

G^(t) = g^(t) 'hj-t^)“® , 

k = lim /2 (x_+h_)^a_ (x_,0) = - 2p lim lim/2 (h -x )^G.{x, ) , 

a X 2 -^-h 2 2 2 2yy 2 o 11 11 

l+2a^(l-6) l-2a^(l-e) 

^o - sinir3 ^ ^ ^ ' 

2 costt 6 + 402(6-!)^ - " ° ' 


“1 = <V2-Vl’/<V^2'^1> 


“2 = (‘^2-^>/<V"l^2’ 


(10) 


where 3 is the power of stress singularity which is real for all material 
combinations, and k is the "stress intensity factor". 

CKdck CKoAAing the InteA^ace. 

In the integral equations (5), If a<hj^, b<h 2 , and 0<c<b, it was shown 
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that the kernels k^^ are bounded in their respective closea- domains [ 2 ]. Re- 
ferring to Figure la and lb, if we let a = hj^, ^”^2' 

the problem becomes one of a crack crossing the interface for which equations 
(5) are still valid. However, in this case following the procedure outlined 
in [ 2 ] and [sl, it can be shown that as Xj^, (i=l,2) and t go to the interface 
together the kernels k^^ (i,j = l,2) become unbounded. The singular parts, 
k^jg(x^,t), of these kernels may be separated by examining the asymptotic be- 
havior of the integrals given by (6). The kernels k. . are given in the Appen- 

1]S 

dix. Together with the kernels (t-x. k. . constitute a set of generalized 

1 xjs 

Cauchy keriic-ls. In this problem the singular behavior of the solution can be 
examined by ejqpressing 

Gl(t) = gj^(t) (hj-t^)"'^ , Q^(t) = g2(t) <h2-t)"^(t-c)“’^ , 0<Re(n,Y)<l . dD 

where g^ and g 2 are again bounded. The characteristic equations to determine 
the constants n and Y obtained from the dominant part of (5) by using the 

complex function technique [3-6] as follows: 

cotirr) = 0 , 

[2cosTry + 4a2(l-Y) ^-aj^-oi2][2cos7rY + 4a^(l-Y)^ 

- [2Y(aj^-a2> + 2-a^ + 3cx2l[2Y(a2-a^) 

+ 2-a^ + 3Ct^] = 0 . (12) 

where the constant (a^, i = l,..,4) are defined in the Appendix. From (12) it 
is found that ri = 0.5, which is the well-known result, and for all material 
combinations there is only one root y satisfying 0<Re(Y)*^l which is always 
real. The same analysis also gives the following relationship between the end 
values of the two bounded functions and g^: 

(h 2 -c)** 2cos7ry + 4ot2 (1-Y) ^- 01^^-012 
^2 ^^2^ (2hi)Y 2Y(a^-a2) + 2 - a^+302 ’ ^l^^l^ 

Equations (12) and (13) are identical to those obtained for the semi-infinite 
planes in [s]. 

A unique solution of (5) in this case too requires two additional condi- 
tions. One condition is given by the single-valuedness of the crack surface 
displacements which, referring to Figure lb and the definitions (3) and (4), 
may be expressed as 

-hi ^^2 
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The second condition is provided by (13) . 

In this problem the stress state is singular at the crack tip 
72 = 0 as well as at the point of intersection of the crack and the interface, 
= + hj^, By examining the stresses around these points, from the 

solution of the problem it can be sho'^m that the stress intensity factors are 
related to the bounded functions g^^ and g 2 as follows [s] . 




2yy' 2' 
U, P 


'*^2 /2 


0 * (h'l ) 1 I 

’'xx “ J ^ '^ixx^^l'y’ “ sin JIX^" (2hi7^ *‘V^ 


g2(h^) 


(h. 


2 2^ r 1-2Y + __1 

2-0) »^1^'^1^2 
V ^ 1 I ^ 1 *.OV 1 

^XV "" T?Y ' ^' u I)C 'u ~ “ iT Tx 'T j" ' ^ 

^2 jV r 1-_2Y _ X ,_ 1 } 


(15) 


In (15) the stress intensity factors k and k are defined explicitly be- 

3CX 

cause of their importance in the study of a possible delamination fracture 
at the interface. However, it should be noted that, because of (13) , in k^^ 
and k^y there is only one independent load factor, say ^2 ^^2^’ 

The solution of the system of singular integral equations with the generalized 
Cauchy kernels subject to the conditions (13) and (14) are obtained by using 
the technique outlined in [ 4 , 6 ]. 

St/LZ6i-^fLZt End PA.obi(>/n 


The stress- free end problem is described in Figure Ic. Analytically the 
problem is simply the limiting case of the crack problem shown in Figure lb 
in which cK). In this case too first the problem for the whole plane (-<»<y<a') 
without the stress-free boundary y * 0 is solved and the homogeneous stresses 
given by (1) are determined. Then a perturbation problem with the end trac- 
tions 

O2yy(x,0) = -Pj^(x) , a2yy{x,0) =■ -P2(x, , 


aixy(x,0) = 0 , a 2 ^y(x, 0 ) = 0 

is considered. The end tractions are statically self-equilibrating and 


( 16 ) 


( 17 ) 


satisfy (see (3)) 

rhi fh2 

j + I p2(x2^‘^*2 ^ ° * 

\ 'hj 

As in most cases, if and P 2 are co* >tant, then 

Pjhj * pjhj - 0 . 


(IS) 


In this problem, basically the integral equations (5) are still valid 
with a = h^, b»h 2 # and c = 0. They may be solved by appropriately separating 
the generalised Cauchy kernels, defining 

and using a technique described in [6], Here, Y is obtained from the charac- 
teristic equation given by (12) and the end points of g^^ and g 2 are related by 


h2 Y 2cos1TY + 4a2(l-Y) -a^-CX2 
^2 ^^2^ * “^h^^ 2Y(aj^-ot2) +T-aj^+3a2 • 


( 20 ) 


Again, to study a possible initiation of delamination fracture, it is useful 
to define and calculate the following stress intensity factors: 


° ii;? , 

y^O sin 
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(2hi) 


J: 2 

• I J 


^2 ^*^2^ r 1-2Y 

1-1 I 4 


>^2-^'^2'^l (2h2)^ ^1^V2 


r_JLi2x_ 

'■Wi+<iy2 


■]> , 


’'xy “ y\xy^^l'y> 


*^1^2 ^l^^l^ j- 1-2Y 


y>0 


cos 


TTY 


P2+*c2yi 


.] _ ^2 ^^2^ |- 1-2Y 


(2h2)'' Wl 


(2h^)^W2 


]} 


^^V2 


( 21 ) 


Nimc/ucaZ Rtiutti 


The material constants used in the nximerical examples and the corresponding 
powers of the stress singularity, 3 for a crack terminating at the interface as 
defined by (10) and Y for a crack crossing the interface (or for the stress- 
free end problem) as defined by (12) and (15) are given in Table 1. Here the 
material combinations A and B are the same and are assumed to appioximate 
boron-epoxy sheets with buffer strips of the same material but different 


stiffness. The material pairs C anu D correspond to Aluminum and Epoxy. Fig- 
ures 2 and 3 give the stress Intensity factor k, for the material combinations 
A and B, respectively# where it is assumed that only material 1 contains a 
crack. Here, k^ is defined by the first equation of (9) for a<h^ and by the 
second equation of (10) for aB>h^. Further results regarding this problem may 
be found in [2]. 


Table 1. The material constants and powers 
of stress singularity 





^^2 

Plane Stress 

Plane Strain 

$ 

Y 

6 

Y 

A 

6.65 

0.33 

0.45 

0.70148 

0.16650 

0.68856 

0.22509 

B 

0.15 

0.45 

0.33 

0.36210 

0.16650 

0.41539 

0.22509 

C 

23.31 

0.30 

0.35 

0.82503 

0.21948 

0.82562 

0.27369 

c 

0.043 

J 

0.35 

0.30 

0.28873 

0.21948 

0.33795 

0.27369 


The stress intensity factors for the crack crossing the interface in the 
material combination A are given in Figures 4,5, and 6. The stress intensity 
factors k , k and k shown in the figures in normalized form are defined 
by (15) . In the examples considered in Figures 2-6 it is assumed that the ma- 
terial is a thin plate (i.e., plane stress case) and is unconstrained in x 
direction. Hence, the constant tractions p^ and P 2 used in the analysis are 
related by 



Figures 7 and 8 show the results for the stress-free end problem described 
in Figure Ic for material con±iinations A and C, respectively. The figures give 
the thin plate results as well as the results for plane strain, i.e., for the 
layered medium which is “thick" in z direction. The stress intensity factors 
shown in these figures are defined by (21). In this case, stress-free boundary 
requires that the constant tractions p^^ and p^ used in the analysis of the per- 
turbation problem satisfy the equilibriuni condition (18) . The results given 

*The problem for the orthotropic strips which is more representative of the 
composite sheets is currently being studied. 



in PiguroB 7 and 8 are normalisod with respuct to » where t«* (h.-fh.) 

The magnitudes of and are determined from the full plane (~»<y<«) under 
the specified external loads. For exaR^Jle, if the medium is subjected to a 

QQ 00 

constant strain e. c. » c in x direction* then 


®l " h2B2(J-Vj)+hj^E^{l-V^> ' ^2 “ " h^ * 


(23) 


If the medium 

in X direction* then 


0<y<«) is subjected to uniform tension 0 


Ixx 2xx 


a 0 


a,h,(ViK2-v,E^) , 

h2E2a-Vj^)+hj^Ej^(l-V2) ' ^2 h^ ^ 


(24) 


Or, if the teinperaturo of the infinite medium (-«<(x,e)<«, 0<y<«) is changed 


by an amount AT, then 

®1®2''2 ^“2““l^ \ 

Pi ■ hjEj(l-Vj)+h^Ej(l-Vj) ' Pj ' - k; Pi 


(25) 


where and are the coefficients of thermal expansion of the layers 1 and 2 
respectively. Input functions for other loading conditions may be obtained in 
a similar way. Needless to say* tho results given in Figures 7 and 8 and those 
which may be obtained on the basis of the equilibrium conditions (17) or (18) 
are applicable to all these loading conditions. 

VlicM&ion and ConcluJ^ioM 


In Figure 2 for a<hj^ and h^'^’O it may be seen that the stress intensity 
factor )c becomes the value obtained from a homogeneoua plane containing peri- 

d 

odic collinear cracks, and for a=hj^ and h^'^’O, as expected, k^-»«. Also* for 

a<h, * as h-***®, k approaches the value obtained for the strip bonded to two 
- 1 2 a 

half planes, which is a special caso of the problem under consideration [7,8], 
These asymptotic values are shown in the Figures by dash-dot lines. Same re- 
sults may be observed in Figure 3 where the crack is in the less stiff ma- 
terial. Here, for a<h, and h.-*0, k^ must and does approach the same periodic 
crack values shown in Figure 2. However, the stress intensity factors in this 
case are much lower than those given in Figure 2 where the crack is in the 
stiff material. 

In Figure 4, for c-*-0, as expected, k ■+<”. It is seen that k also 

C 

when c'+hj. The reason for this is that for c^h^ the power of stress singu- 
larity 0 is greater than 0.5 (material combination A), whereas k shown in 

c 
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Fi>juru 4 i« calculated on the hauls of one half t»ower singularity. It may also 
bu observed that becomes greater as h^/h^ Increases. Since the power of the 
stress singularity 0 f- r a crack terminating at the interface is always greater 
than Y for a crack crossing the interface# k and k also go to infinity as 

aX 

the crack tip approaches the interface (Figures 5 and 6) . 

When the singular point is at the Interface as in problems for a crack 
touching or crossing the interface# since the power of singularity is not 0.5, 
for a crack propagation initiating at this point the stress state around the 
crack tip does not remain self-similar. Therefore# in this case most of the 
conventional fracture and fatigue models do not seen, to be applicable. In the 
absence of a physically more acceptable criterion# at present one may assume 
that for this type of fracture problems the ''maximum stress criterion" is 
adequate. This criterion may bn stated as "the fracture propagation will taka 
place radially in the direction 0 «0^ for which the cleavage stress Cqq is 
maximum and when 


#0) “ constant# (26) 

yo p - c 

where 6 is the siao of the fracture process zone around the crack tip". The 
P 

process zone size depends on the microstructure and continuum pi’opertics of 

the material as well as on the environmental conditions. The critical stress 

a is considered to be a "material constant" representing the cohesive 
c 

strength of the constituent materials and, for the appropriate 6 , is deter- 
mined from controlled experiments for an idealized geometry and loading. If 
the "weak link" around the singular point is the Interface, then (26) may be 
modified as follows 


((7 ) 
' c 


12 


,2 

06 

.< 3 ^, 1 '’ , 

“ee =■ “ 

IG 

* “rO • 

“00 ‘ ° 


(27) 


where (o^) represents the adhesive strength of the bond, fd^Q is the friction 
resistance# and the scresscs Oyg and are calculated at a distance r= 

from the singular point and represents tne process zone size for the particular 
joint. 

Finally, from Table 1 and Figures 7 and 8 it may be observed that in the 


stress-free end problem the stress intensity factors# k , k and the power 
of the stress singularity Y for the plane stress conditions are greater than 
those for the plane strain cate. This moans that for the stress-free end 


problem shown in Figure Ic the atreaa utato around the bonded cornor« under 
plane atreaa conditions ia expected to bo much more severe then that under 
plane strain conditions. 
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APPENDIX 


Rafarring to (5) and (6), let the asymptotic expression of for s-*<» 


he i»e*# 


Kij(Xi,t,s) - K^^^(Xj^,t,s) ♦ K^^^(x^,t,s) , (i,j-l,2) . 


(Al) 


Then defining 


>Ci^(Xi,t) - ' (i,j-l,2) , 


(A2) 


the singular parts* bounded parts* ^he kernels can be 

obtained by substituting from (Al) and (A2) into (6) . Thus* using the ex- 
pressions for given in Cs]* the singular kernels a*"® found to bet 

kll.lXj.t) . ^ 6aj(h,*»j) ai;-* 

j'Vl’ 5^ ♦ 2“j<h,-xil' ^ 

1 v*jv j 


3d *d 

+ [-4-^ - 6a 


1-1 


2-a +3a _ 

>tl2s^*l'^) - [(a^-a^) (hj^-x^) + 5 ^ 

2-a +3a 

+ [(a^-aj^) (hj^+Xj^) ^ + 5 ][t-(hj+h 2 -*-Xj^) 3 * 

2 -a -«3a _ 

*‘21s^’‘2'^* ” f *^2"*2^ * 2 3[t- (hj^+h^-x^) 3 


’*22s^’‘2'^^ " ^ 2 ~ dx ■*■ ^“4^^2"’‘2^ dx 3[t- (2h2-X2> 3 

* M 


-1 




(A3) 
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Figure 1. The geometry of bonded strips and the stress-free end. 





1 .0 / 0.9, r 0.8, r 0.7, r 0.6, r 0.5, r 0.4, r 0.3, r 0.2 


Figure 2. 


Stress intensity factor in bonded strips for the crack in the 
stiffer n«terial (Ui»6.65U2. Vi.0.33. V2-0.45; the scale w 
the right corresponds to the broken strip case, i.e., to a -hi 
dash-dot lines give the asymptotic values for h 2 -*®) • 






Figure 3. Stress intensity factor in bonded strips for the crack in less 

stiff material (yi*6.65U2. Vj^»0.33, V2*0.45; the scale on the 
right corresponds to b = h 2 > . 











Figure 5. Normal component of the stress intensity factor at the inter 
section of the crack and the interface for the crack crossing 
the interface where material 1 is fully cracked (Wi«6.65U2» 
Vj^»0.33, V 2 -O. 45 , f»hj^+h2-c). 
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Figure 6. Shear component of the stress intensity factor at the inter- 
section of the crack and the interface for the crack crossing the 
interface where material 1 is fully cracked 6.65^2* Vj^*0.33» 

V2*0.45, -d *= h]^+h2~c) . 
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Fiqure 7. Normal and shear components of the stress intensity factor, kj^j^ 
and k^y at the intersection of the interface and the stress-free 
boundary in a layered half planet W], * ^1 * V2*0.4S, 



Normal and shear components of the stress intensity factor, k^x 
and kxy at the intersection of the interface and the stress-free 
boundary in a layered half plane; Vx~0.3, V2*0.35, 

L ■ h2^+h2 • 


Figure 6. 


